An empirical test for cellular automaton models of traffic flow 
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Based on a detailed microscopic test scenario motivated by recent empirical studies of single- 
vehicle data, several cellular automaton models for traffic flow are compared. We find three levels 
of agreement with the empirical data: f ) models that do not reproduce even qualitatively the most 
important empirical observations, 2) models that are on a macroscopic level in reasonable agreement 
with the empirics, and 3) models that reproduce the empirical data on a microscopic level as well. 
Our results are not only relevant for applications, but also shed new light on the relevant interactions 
in traffic flow. 
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I. INTRODUCTION 

For a long time the modeling of traffic flow phenom- 
ena was dominated by two theoretical approaches (for a 
review, see e.g., [1-4]). The first type of models, the so- 
called car-following models, are based on the fact that the 
behavior of a driver is determined by the leading vehicle. 
This assumption leads to dynamical velocity equations 
which in general depend on the distance to the leading 
vehicles and on the velocity difference between the lead- 
ing and the following vehicle. An alternative approach, 
which is also well established in traffic research, does not 
treat the individual cars but describes the dynamics of 
traffic networks in terms of macroscopic variables. Here 
traffic flow phenomena are treated in analogy to the dy- 
namics of compressible viscous fluids. 

Both approaches are still widely used by traffic engi- 
neers, but for practical purposes they are often not suit- 
able. One of the main problems of present car-following 
models (e.g., see [5-9]) is that they are difficult to treat 
in computer simulations of large networks. On the other 
hand also the macroscopic approaches lead to some diffi- 
culties although large networks can be treated in prin- 
ciple. First of all, present macroscopic models use a 
large number of parameters which have partly no coun- 
terpart within empirical investigations. In addition to 
that, the information that can be obtained using macro- 
scopic models is incomplete in the sense that it is not 
possible to trace individual cars. 

In order to fill this gap cellular automaton (CA) models 
have been invented [10,11]. CA models are microscopic 
models which are by design well suited for large-scale 
computer simulations. A comparison of the simulations 
with empirical data shows that already very simple ap- 
proaches give meaningful results. In particular they can 
be used in order to simulate dense networks like cities [12] 
which are controlled by the dynamics at the intersections. 
For highway traffic, however, a more detailed description 
of the dynamics seems to be necessary. 

In this work we want to discuss the realism and the 
limitations of a number of CA models. Our choice is re- 
stricted to models that are discrete in space and time, 



which e.g. excludes the approach by Krauss et al. [13], 
and have local interactions only, excluding models as 
the Galilei- invariant model introduced in [14]. We com- 
pare simulations of the CA model proposed by Nagel 
and Schreckenberg, that is to date the most frequently 
used CA approach for traffic flow, the VDR model [15] 
which realizes a so-called slow-to-start rule, the TOCA- 
model of Brilon et al. [16], the model of Emmerich and 
Rank [17] based on the use of velocity-gap matrices, and 
the approach by Helbing and Schreckenberg [18] which 
represents a model with a more sophisticated distance 
rule. Finally we discuss the recently introduced brake 
light model [19, 20] that was suggested in order to give a 
reliable reproduction of the microscopic empirics and the 
model by Kerner, Klenov and Wolf [21], focusing more on 
the macroscopic properties of the three phases of traffic 
flow. 

We will compare the ability of these models to repro- 
duce the empirical findings. This requires using a mea- 
surement procedure in the simulations which models the 
detectors on the highway. Analogous to the empirical 
setup of [22] the simulation data are evaluated by a vir- 
tual inductive loop, i.e., speed and time-headway of the 
vehicles are measured at a given link of the lattice. The 
measurement process is applied after the update of the 
velocity has been carried out, but right before the move- 
ment of the vehicles. This implies that the gap to the 
preceding vehicle does not change significantly during 
the measurement. These simulation data are analyzed 
regarding individual and aggregated quantities, as it has 
been done in recent empirical investigations [22-24] . 

Although most of the empirical data sets have been 
collected at multi-lane highways, we have performed our 
simulations on a single-lane road in order to reduce the 
number of adjustable parameters. This approach is jus- 
tified because the empirical data sets are selected such 
that multi-lane effects are of minor importance. They 
might play a role for synchronized traffic of type (i) and 
(m), as it has been recently argued in [25], but in any 
case these types of synchronized traffic are much rarely 
observed than synchronized traffic of type (Hi) [22, 24]. 

We will also not consider effects by a mixture of diffcr- 
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ent vehicle types, e.g. there are no trucks in our simula- 
tions. The fraction of slow cars has not been determined 
from the empirical data. Furthermore these data have 
been collected on a highway with speed limit such that 
disorder effects through slower cars are expected to play 
a minor role. We believe that inclusion of disorder will 
not change our results qualitatively, but can lead to a 
better quantitative agreement in some cases. 

Before we start the analysis of the above mentioned 
CA models, we will introduce an empirical test scenario. 
It will be microscopic and local to make it easily compa- 
rable to online data provided e.g. by inductive loops. In 
contrast, the detection of complex spatio-temporal struc- 
tures [26] is more difficult to achieve in an automated 
way. It would require the investigation of interface dy- 
namics whereas in our scenario only bulk properties are 
studied. This test scenario also verifies the reproduc- 
tion of empirical traffic states on a microscopic level, a 
task that cannot be fulfilled by macroscopic models. The 
empirical results have been chosen with respect to their 
reproducibility and the ability to distinguish between the 
different states of traffic. This scenario will be discussed 
in the next section. 



II. EMPIRICAL FACTS 

In order to probe the accuracy and the degree of re- 
alism of the different models one has to introduce a test 
scenario that includes the most important empirical find- 
ings. The difficulty in defining such a scenario is due to 
the fact that the empirical results may depend strongly 
on the particular environment. Therefore one has to try 
to extract the results that really characterize the behav- 
ior of the vehicles. As an additional difficulty mostly 
aggregated data have been analyzed which are known to 
be largely dependent on the road conditions, e.g., the ca- 
pacity of an upstream bottleneck. A number of results, 
however, is of general nature as we will discuss below. 

Even more conclusive are empirical investigations that 
use single- vehicle data. These measurements can be com- 
pared directly to the simulation results and include im- 
portant information concerning the microscopic structure 
of vehicular traffic. Unfortunately only a small number 
of empirical investigations based on single-vehicle data 
exists so far. Our discussion refers to the empirical stud- 
ies of refs. [22-24] . In particular, in order to reduce the 
effects of disorder, the results of [22] (except for the time- 
headway distributions, see below) are used for the com- 
parison with simulation data. These data have been col- 
lected on a highway where a speed limit applies. This 
facilitates the comparison with modeling approaches. 

The empirical findings that are taken as a basis for 
the comparison with the model results have been ob- 
tained from inductive loops. Measurements by induc- 
tive loops, which represent the most frequently used mea- 
surement devices, give information about the number of 



cars passing, their velocities and the occupation times. 
These direct measurements are also used in order to cal- 
culate other quantities, e.g., the spatial distance d n via 
d n = v n -\th (where v n -\ denotes the velocity of the pre- 
ceding car n — l,th the time- headway between car n — 1 
and car n). 



A. Temporally aggregated data 

The most important empirical quantity is the relation 
between the averaged observables flow and density, i.e., 
the fundamental diagram. There exists a longstanding 
controversy (see e.g. [26, 27] and references therein) about 
the "correct" functional form of the fundamental diagram 
and a large number of possible forms have been suggested 
to be compatible with empirical data [28]. A more con- 
sistent picture was established after the work of Kerner 
and coworkers who distinguished at least three different 
phases of traffic flow [29], i.e., free flow, synchronized traf- 
fic and wide jams, that have to be analyzed separately. 
We will follow this scheme and summarize the empirical 
findings accordingly. 

Usually these measurements are stored as averaged val- 
ues of certain time-intervals. We discuss results for the 
fundamental diagram, i.e., the flow density relation, in 
the different traffic phases that are based on one-minute 
data. The results for the functional form of the flow are 
shown, as far as possible, in dependence of the spatial 
density p(t). The density can be calculated from 
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where J(t) denotes the number of cars passing the de- 
tector with an average velocity v(t) in the corresponding 
time interval. 

Free flow traffic is characterized by a large value of the 
average speed. One basically observes two qualitatively 
different functional forms of the fundamental diagram, 
i.e., that the linear regime extends up to the observed 
maximum of the flow or that one has a finite curvature 
in particular for densities slightly below the density of 
maximum flow [24, 30] . The finite curvature is a conse- 
quence of an alignment of speeds, i.e., close to the optimal 
flow it is not longer possible to drive systematically faster 
than the trucks. This point of view is supported by the 
empirical results taken from highways where a quite re- 
strictive speed limit is applied that can be reached even 
by trucks [22]. 

In this case the whole free flow branch is linear. For 
our purposes the linear form of the fundamental diagram 
is relevant, because we use a single type of cars in the 
simulations, with a maximal velocity that is given by the 
slope of the free flow branch. When simulating a section 
of the highway where no speed limit is applied, one has to 
take a distribution of maximal speeds. This distribution 
can be obtained from the empirical velocity distributions 
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at very low densities, where interactions between cars can 
be neglected. 

In the congested regime one distinguishes between syn- 
chronized traffic and wide jams. In the synchronized 
phase, the mean velocity of the vehicles is reduced, com- 
pared to the free flow, but the flow can take on values 
close to the maximum flow. Moreover, strong correla- 
tions between the density on different lanes exist caused 
by lane changings. 
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FIG. 1. Time-traced fundamental diagram of the two con- 
gested states (from [30]). Synchronized traffic is characterized 
by strong fluctuations of the density and flow. The measure- 
ments for wide jams are similar to measurements in free flow 
but with much smaller average velocity. 



The synchronized state has been subdivided into three 
types, which differ in the characteristics of the time se- 
ries of density and flow: In synchronized traffic of type 
(i) constant values of the density and the flow can be 
observed during a long period of time. In synchronized 
traffic of type (ii) the flow depends linearly on the den- 
sity similarly to free flow, but the mean velocity is re- 
duced considerably. In synchronized traffic of type (Hi) 
irregular patterns of flow and density can be observed 
(Fig. I). In our article we concentrate on synchronized 
flow of type (Hi), because the two other types of syn- 
chronized traffic have been rarely observed and it is not 
confirmed whether they are generic phases of traffic flow. 
An identification of synchronized traffic by means of the 
fundamental diagram may be misleading, because the re- 
sults often depend on the averaging procedure. A more 
sensitive check is to identify the different types of traffic 
states by means of the cross-correlation cc(p, J) of the 
density p and the flow J [22]: 



with A A = (A 2 ) — (A) 2 denoting the variance of the 
observable A. The linear dependency of the flow and the 
density in the free flow state as well as in the wide jam 
state leads to cross-correlations of « 1, whereas irregular 
patterns of the flow and the density in the synchronized 
traffic of type (in) lead to cross-correlations of w 0. 

It is worth pointing out that the notion of "synchro- 
nized traffic" is still very controversial [26,27]. We em- 
phasize here, that we use an objective criterion, namely 
the vanishing of the cross-correlation function (2) for the 
classification. Within the empirical single-vehicle data 
sets available the other two synchronized states could 
not be clearly identified. Therefore it was not reason- 
able to include these states into the test scenario. The 
characteristics of synchronized traffic of type (Hi), how- 
ever, have been clearly distinguished from free flow and 
jammed states by the criterion cc(p, J) w 0. Therefore 
any detailed model should be able to reproduce this class 
of synchronized states. 

Fig. 1 includes a typical measurement of the fundamen- 
tal diagram that correspond to wide jams. Surprisingly 
these measurements reveal quite small values of the den- 
sity, although the road is almost completely covered by 
cars. This seemingly incorrect result is due to the local 
nature of the measurement (see [22] for a detailed dis- 
cussion). Thus, the form of the fundamental diagram in 
the jammed state is similar to free flow traffic, but with 
a small average velocity. 

The jammed branch of the fundamental diagram is of- 
ten not reproduced by CA models, because they use the 
inverse density of a jam in order to calibrate the unit of 
length. Within these approaches jams are compact. In 
this case (almost) no internal flux is observed. The mod- 
eling of jams can however be meaningful, if the upstream 
velocity and other macroscopic characteristics of jam are 
reproduced. 



B. Single vehicle data 

Nowadays some empirical studies exist that have an- 
alyzed single- vehicle data from counting loops [22-24]. 
These studies are of great importance for the modeling of 
traffic flow because they give direct information about the 
"microscopic structure" of traffic streams. The data usu- 
ally include direct measurements of the time-headways 
and the velocities of the vehicles as well as the occupa- 
tion time of the detector. Similar to the time-averaged 
observables the results for the microscopic quantities dif- 
fer qualitatively in the different phases. 

The first quantity we look at is the time-headway dis- 
tribution *, i.e., the time elapsing between two cars pass- 



cc(p, J) 



(p(t)J(t + r))-(p(t))(J(t + r)) 
^Apjt) ^AJ(t + r) 



(2) 



* Since the time-headway distribution of [22] in free flow as 
well as in the synchronized state shows some peculiarities due 
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ing the detector. This quantity is the microscopic ana- 
logue to the inverse flow. In free flow traffic one has found 
that the distribution at short times and also the position 
of the maximum is independent of the density (Fig. 2). 

The cut-off at small time-headways as well as the typ- 
ical time-headway in free flow traffic are important ob- 
servables which have to be reproduced by the microscopic 
models. The exact shape of the distribution may also de- 
pend on the relative frequency of slow vehicles, because 
this determines the fraction of interacting vehicles at a 
given density. 

The time-headway distributions in synchronized traf- 
fic ' differ systematically from the free flow distributions 
(Fig. 3). In synchronized traffic the distributions have 
a maximum that is much broader than that in free flow 
traffic. The maximum is less pronounced and its position 
depends significantly on the density. 




FIG. 2. Empirical time-headway distributions, i.e., the rel- 
ative frequency of a given time-headway, in free flow traffic. 
The distributions are normalized, i.e., ^ P(th)-Ath = 1. The 
data are classified in different density regimes by the corre- 
sponding one-minute data of the density. For a given road 
section one obtains a maximum that is independent of the 
density and a minimal headway of 0.2 s. 



to an error of the measurement software [24] , new measure- 
ments at the same location have been conducted. 

t Unfortunately, new measurements taken from the detector 
location used in [22] do not provide a sufficient amount of 
data of the synchronized state. Since the time-headway dis- 
tributions of [22] cannot be used, the distribution is calculated 
from data sets taken from [24]. This is justified because the 
effects of a speed limit can be neglected at larger densities. 
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FIG. 3. Same as Fig. 2 but for synchronized traffic. The 
functional behavior of the distribution at short times depends 
on the density. 

In the presence of wide jams one has to distinguish be- 
tween the jam itself and its outflow region. In the jam one 
finds evidently a broad distribution of time-headways, 
because cars are blocked for quite long times. In the 
outflow region of a jam, however, one observes that the 
typical time-headway is of the order of ~ 2 s. 

The characteristics of traffic jams are one of the ex- 
tensively studied phenomena in traffic flow. Wide traffic 
jams can be identified by a sharp drop of the velocity and 
the flow to negligible values in the time-series. Moreover 
traffic jams move upstream with a surprisingly constant 
velocity (typically 15 km/h [31]). The upstream velocity 
is intimately related to the outflow J out from a jam which 
also takes on constant values for a given situation. This 
allows the observed coexistence of jams. The coexistence 
is facilitated because the outflow from a jam is consider- 
ably smaller than the maximal flow J max , such that no 
new jams emerge in the outflow region of a jam. Empir- 
ically one observes the ratio J max / </ ut ~ 1-5 [32]. The 
outflow and the upstream velocity of a jam can therefore 
also be used to calibrate the model. The precise data for 
the average upstream velocities and J out may also serve 
to evaluate the average space / that is occupied by a car 
in a jam. Usually I, and not the average length of the 
vehicles, represents the length of a cell in the CA mod- 
els. I may also be used to assign a reasonable value of 
the velocity of cars in a jam, i.e., v n = l/th- 
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FIG. 4. Empirical optimal-velocity (OV) functions, i.e., 
speed-distance relations. The figure shows the mean velocity 
for a given spatial distance in free flow and congested traffic 
for different densities. 



The final test of the models comes from the velocity 
distance relation in the different traffic phases (Fig. 4). 
This relation, also called optimal- velocity (OV) function, 
characterizes in great detail the microscopic structure of 
the different phases. Some models use OV-curves directly 
as an input [33]. In any case this quantity is a sensitive 
test concerning the reproduction of the microscopic struc- 
ture of highway traffic. In the free flow regime the asymp- 
totic velocity does not depend on the density, but is given 
by the applied speed limit. In the congested regime this 
asymptotic velocity is much smaller than in free flow, i.e., 
cars are driving slower than the distance-headway allows. 
This is a direct effect of the vehicle-vehicle interactions 
[22] and should therefore be reproduced by any realistic 
traffic-model. 



III. SIMPLE STOCHASTIC CA MODELS 

Throughout this article we investigate microscopic 
traffic models that are discrete in space and time. The 
discreteness of the model has the advantage of allowing 
direct and very efficient computer simulations, and in 
particular without any further discretization errors. The 
discreteness of the model, however, also leads to some 
difficulties, in particular when describing congested traf- 
fic. E.g., in congested traffic a continuous range of typi- 
cal velocities exists that depend strongly on the density. 
This velocity interval is mapped on a discrete set of ve- 
locity variables. So even for an optimal reproduction of 
the traffic state an upper limit for the accuracy of the 
model exists. Therefore, one has to find a compromise 



between the degree of realism and the level of complexity 
by choosing an appropriate discretization of the velocity. 

Moreover, the temporal discretization introduces a 
characteristic time-scale. This time-scale can be under- 
stood, if a parallel update is applied, as the effective 
reaction-time of the drivers, which is included explicitly 
in car-following models. Furthermore, the temporal dis- 
cretization becomes obvious as peaks in the measurement 
of the time-headways th- The finer the discretization the 
less pronounced the peaks. In order to increase the res- 
olution the time-headways in the simulations are calcu- 
lated via the relation th = - with the velocity v of the 
vehicle and the distance-headway d to the preceding vehi- 
cle. Nevertheless, the minimal resolution is restricted by 
the discretization that determines the minimal th differ- 
ence in free flow — - — with the length I and the maximum 
velocity w max of a vehicle. In order to facilitate a com- 
parison with the empirical time-headway distribution the 
distributions are normalized via ^P{th) ■ At/, = 1- 

Below we discuss a number of traffic models in detail 
and with respect to their agreement with the empirical 
findings of our test scenario. Beyond that we demand 
that each model reproduces some basic phenomena, like 
spontaneous jam formation, and fulfills minimal condi- 
tions as, e.g., being free of collisions. These conditions 
are generally understood as fulfilled, if the opposite is 
not explicitly stated. In particular, deterministic models 
(e.g. [34,35]) are not a subject of this study. They can 
not reproduce the spontaneous formation of jams [11] 
which are the result of an inherent stochasticity of traffic 
flow rather than a consequence of perturbations. 

Our simulations are performed on a periodic single- 
lane system. This simple structure of the system is in 
sharp contrast with realistic highway networks. It is nev- 
ertheless justified, because it has been shown for a large 
class of models that different boundary conditions se- 
lect different steady states rather than changing their mi- 
croscopic structure [36]. Therefore the boundary condi- 
tions are of great importance if one tries to reproduce the 
spatio-temporal structure on a macroscopic level. How- 
ever, in comparison with local measurements an appro- 
priate traffic model should be able to reproduce the em- 
pirical results also if periodic boundary conditions arc 
applied. Furthermore the restriction to a single lane is of 
minor importance for the empirical test scenario which 
has been discussed in the previous section. In the sim- 
ulations system sizes of L > 10000 cells have been used 
which is sufficient to reduce finite-size effects. Typical 
runs used 50000 time steps to reach the stationary state 
and measurements. 

We also want to emphasize that for each model all 
simulations have been performed with a single set of pa- 
rameters. Some of the model parameters can be directly 
related to a given empirical quantity. In this case we have 
chosen the value that leads to an optimal agreement with 
the related observable to avoid ranking the importance of 
the empirical findings. For a particular application of the 
model, however, the reproduction of a certain quantity 
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might be of special interest and therefore a calibration of 
the model different from ours might be more appropriate. 



A. The CA model of Nagel and Schreckenberg 

The model introduced by Nagel and Schreckenberg [10] 
(hereafter cited as "NaSch model" ) is the prototype of 
microscopic models that we discuss. The important role 
of this model is mainly due to its simplicity which allows 
for very fast implementations. In fact the NaSch model 
is a minimal model in the sense that every further sim- 
plification leads to a loss of realism. We will also use it 
as a reference for other models that will be introduced 
by giving the relation to the NaSch model. 

The NaSch model is a discrete model for traffic flow. 
The road is divided into cells that can be either empty or 
occupied by car n with a velocity v n = 0, 1, n max . Cars 
move from the left to the right on a lane with periodic 
boundary conditions and the system update is performed 
in parallel. 

For completeness we repeat the definition of the model 
that is given by the four following rules (t < t\ < t 2 < 
t+l): 

1. Acceleration: v n (ti) = min{w„(t) + l,w max } 

2. Deceleration: v n {t2) = min{w„(ti), d n (t)} 

3. Randomization: v n (t+l) = max{i>„(t2)~ 1, 0} with 
probability pdcc (otherwise v n (t + 1) = v n {t,2)) 

4. Motion: x n (t + 1) = x n (t) + v n (t + 1) 

with the velocity v n , the maximum velocity w max and the 
position x n of car n. d n (t) specifies the number of empty 
cells in front of car n at time t. 

For a given discretization the model can be tuned sim- 
ply by varying the two parameters v maK and pdcc- The 
value of "Umax mainly affects the slope of the fundamen- 
tal diagram in the free flow regime while the behavior in 
the congested regime is controlled by the braking noise 
Pdoc- Each time-step At corresponds to 1.2 s in real- 
ity in order to reproduce the empirical jam velocity at 
a given cell length of 7.5 m. The length of a cell corre- 
sponds to the average space occupied by a vehicle in a 
jam, i.e., its length and the distance to the next vehicle 
ahead. This choice is in accordance with measurements 
at German highways on the left and middle lane, where 
the density of trucks is low [32]. Due to the parallel up- 
date an implicit reaction time is introduced which has 
to be considered when choosing the unit of time. This 
time is not the reaction time of the driver (that would 
be much shorter) but the time between the stimulus and 
the actual reaction of the vehicle. The value we have 
chosen allows to reproduce the typical upstream velocity 
of a jam. 




We tune the two free parameters of the model by ad- 
justing the slope in the free flow regime and the maxi- 
mum of the fundamental diagram. Fig. 5 shows the re- 
sulting fundamental diagram using w max = 112 km/h = 
5 cells/timestep and pdec = 0.16 which has to be com- 
pared with the empirical results. 
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FIG. 5. Local fundamental diagram of the NaSch model for 
fmax = 112 km/h = 5 cells/At, At — 1.2 s and p doc = 0.16. 
A cell has a length of 7.5 m. 



By tuning the parameters we could reproduce quite 
well the free flow branch of the fundamental diagram: 
Both, the slope as well as the maximum is in agreement 
with the empirical findings. For congested traffic, how- 
ever, the model fails to reproduce the two distinct phases, 
in particular the characteristics of synchronized traffic 
are not matched. This interpretation of the flow data is 
supported by measurements of the cross-correlation func- 
tion that is negative in the corresponding density regime. 
In the presence of wide jams the flow is proportional to 
the densities as found by empirical observation. But also 
for wide jams differences exist. In real measurements the 
branch extends up to quite large densities (~ 70 veh/km), 
while the simulation results are restricted to lower den- 
sities (~ 40 veh/km). 

Next we discuss the model on a microscopic level. As 
mentioned above the upstream velocity of wide jams can 
be tuned by choosing the appropriate discretization At of 
the time. We have verified our calibration by initializing 
the system by a large jam and measuring the velocity of 
the upstream propagation of the jam front. As expected 
our result is in agreement with the empirical data. Nev- 
ertheless the dynamics of jams in the NaSch model is in 
contradiction to empirical findings since its outflow from 
a jam equals the maximal possible flow. This implies 
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that the observed parallel propagation of jams cannot be 
reproduced by the NaSch model. 

The time-headway distributions of the NaSch model 
(see also [37]) also mismatch with empirical data (Fig. 6). 
Due to the discreteness of the model and the unique max- 
imal velocity of the cars the distribution function has a 
peaked structure^. 
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FIG. 6. Normalized time-headway distribution of the 
NaSch model in free flow and congested traffic for different 
densities. 



But more important than that is the absence of time- 
headways shorter than the chosen unit of time. This 
implies that we cannot reproduce the cut-off at short 
times and the upstream velocity of jams at the same time. 



*The time-headway distributions have a resolution that is 
finer than the unit of time which was assigned to an update 
step. This is possible because we calculate the exact passing 
time of the car from its position and velocity after executing 
the time-step. An example of a direct measurement of time- 
headways can be found in [37]. 
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FIG. 7. OV function in free flow and congested traffic of 
the NaSch model for different densities. 



Finally, we also discuss the optimal velocity curves of 
the model (Fig. 7). In congested traffic one observes only 
a very weak dependence of the "optimal velocity" on the 
density. This is due to the short range of interactions 
in the model and the strong acceleration of the cars. So 
we neither observe a significant density dependence nor a 
sensitivity to the traffic state. This is a serious contradic- 
tion to the empirical findings, related to an incomplete 
description of the microscopic structure of the model. 



B. VDR model 

A step towards a more realistic CA model of traf- 
fic flow was done by the so-called velocity-depcndcnt- 
randomization (VDR) model [15] that extends slightly 
the set of update rules of the NaSch model. In this 
model, a velocity-dependent randomization Pd ec (v) is in- 
troduced that is calculated before application of step 1 
of the NaSch model. As simplest version, a different Pd ec 
for cars with v = was studied: 



Pdcc(w) 



P0 



for v = 



p for v > 



(3) 



with po > P (slow-to-start rule). 

The additional rule of the VDR model has been intro- 
duced in order to reproduce hysteresis effects. This is 
indeed possible, because the new parameter p allows to 
tune the velocity and outflow of wide jams separately. As 
a side effect it is now possible to reproduce the observed 
short time-headways by keeping the unit of time small 
and the empirical observed downstream velocity of jams. 
The parameters of the model were chosen in the following 
way: The unit of time was adjusted in order to match the 
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position of the maximum of the time-headway distribu- 
tion. Then we have chosen the parameter po such that we 
could reproduce the measurements of the upstream ve- 
locity of a jam. Finally the values of w max and p ensure 
a good agreement in the free flow branch. The behav- 
ior found in the VDR model is typical for models with 
slow-to-start rules [38,39]. 
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quite large positive values of the cross-correlation cc(J, p) 
of the density and the flow. 
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FIG. 9. Normalized time-headway distribution of the VDR 
model in free flow and congested traffic for different densities. 
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FIG. 8. Local fundamental diagram of the VDR model for 
Umax = 108 km/h = 3 cells/At, At = 0.75 s, p = 0.58 and 
p = 0.16. 

Fig. 8 shows the local fundamental diagram of the 
VDR model. For the parameter values obtained by the 
above procedure only very weak hysteresis effects are ob- 
served. Obviously the model fails to reproduce the em- 
pirically observed congested phase correctly. Compared 
to the NaSch model the mismatch of the fundamental di- 
agram in the congested regime is even more serious, i.e., 
we cannot identify at all a density regime as synchro- 
nized traffic. The reason for this is a stronger separation 
between free flow and wide jams, which are compact. 
Therefore one does not observe any flow within a jam if 
a stationary state of a periodic system is analyzed. In 
case of open boundary conditions a slight broadening of 
the free-flow branch has been observed, if the detector 
is located close to the exit of the highway section. This 
effect is due to the smaller length scale of jams close to 
the exit, which leads to a larger weight of accelerating 
cars. Due to the coarsening of the jam size this effect 
vanishes in the bulk of the system [40,41]. 

But in any case, this way of generating synchronized 
states by the boundary conditions does not agree with 
the empirical situation, because one cannot reproduce the 
large spatial and temporal extension of the synchronized 
state. The missing synchronized traffic phase leads to 
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FIG. 10. OV function of the VDR model in free flow and 
congested traffic for different densities. 

The time-headway distribution of the VDR model 
differs in two points from the empirical observations 
(Fig. 9). (i) The unit of time is a sharp cut-off, i.e., 
the short time characteristics of the time-headway distri- 
bution is not in agreement with the empirical findings, 
(ii) We do not observe a density dependence of the maxi- 
mum in congested traffic. Similar results are obtained for 
the OV functions, that do not depend on the density or 
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the traffic state (Fig. 10). This result is a consequence of 
the microscopic structure of high density states. At large 
densities compact wide jams and zones of free flow traffic 
coexist, separated by a narrow transition layer. Now, our 
virtual "detector" measures only moving cars and there- 
fore almost freely moving cars even at large densities. 

The major achievement of the VDR model is the cor- 
rect description of the dynamics of wide jams which is 
similar to the so-called local cluster effect [42] found in 
hydrodynamical models. The outflow from a jam is lower 
than the maximal flow, and therefore jams do not emerge 
in the outflow region. This effect leads to the increased 
stability of jams, including the empirically observed par- 
allel upstream motion of two jams. 

The analysis of the VDR model showed even more 
clearly the effect of a missing synchronized traffic phase. 
While in the NaSch model the density can be chosen 
such that a scattered structure in the fundamental dia- 
gram appears, we obtain rather pure free flow states and 
wide jams for the VDR model. Contrary the VDR model 
gives a much better description of the dynamics of jams. 
In contrast to the NaSch model, the VDR model is able 
to reproduce e.g. the parallel motion of coexisting jams 
[40,41]. Although this phenomenon is rarely observed it 
should be reproduced by a realistic traffic model, because 
it is a sensitive for the correct description of the motion 
of jams. In case of the NaSch model this pattern is not 
observed, because new jams can form in the downstream 
direction of a jam. 



C. The time-oriented CA model 

Based on the CA model of Nagel and Schreckenberg, 
Brilon et al. [16] proposed a time-oriented CA model 
(hereafter cited as TOCA) that increases the interaction 
horizon of the NaSch model (where cars interact only for 
d < v) and therefore changes the car-following behavior. 

Compared to the NaSch model the acceleration step is 
modified, i.e., a car accelerates only if its temporal head- 
way th = d(t)/v(t) is larger than some safe time-headway 
t s . But even for sufficiently large headways the accelera- 
tion of a vehicle is not deterministic, but is applied with 
probability p ac . As a second modification also the ran- 
domization step is modified, i.e., it is performed only for 
cars moving with short time- headways (th < t s ). The 
limited interaction radius of this third step leads, for a 
given value of pdec, to a reduction of the spontaneous jam 
formation. 

The update rules then read as follows (t < t\ < t 2 < 
t+l): 



1. if (t h > t s ) then 

v n (ti) = min{v n (t) + l,v„ 



} with probability p a 



2. v n (t 2 ) 



i{v n (ti),d n (t)} 



3. if (th < t s ) then 

v n (t+l) = max{w„(i 2 ) — 1, 0} with probability pde 



4. X n (t + l) = x n (t) + v n (t + 1) 

with t s = 1.2, p ac = 0.9 and pdec = 0.9 [16]. For the com- 
parison with the NaSch and VDR model we use w max = 4. 

With this choice of t s the update rules can be simplified 
for w max < 4, because of the discrete nature of the model: 

1. v n (ti) = mm{v n (t) + l,v max } with probability p ac 

2. v n (t 2 ) = mm{v n (ti) , d n (t)} 

3. if (v n (t + l) <d n (t)) 

v n (t + 1) = max{d n (t) — 1, 0} with probability pdec 

4. X n (t + l) = x n (t) + v n (t + 1) 
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FIG. 11. Fundamental diagram of the TOCA model. As 
discretization we used a cell length of 7.5 m and a time-step 
corresponding to At = 1 s in reality. The parameters of 
the model are chosen as t s = 1.2, p ac = pdoc = 0.9 and 
«max = 4 cells/At = 108 km/h. 



As expected for this parameterization of the model we 
obtain results for the fundamental diagram that are sim- 
ilar to the NaSch model (Fig. 11). 

The absence of spontaneous velocity fluctuations at 
low densities, however, implies that pdec has to be chosen 
quite large in order to obtain realistic values of the max- 
imal flow. At the same time large values of the braking 
probability lead to the formation of jams at low densi- 
ties, such that it is difficult to obtain density fluctuations 
with amplitudes comparable to the empirically observed 
values. 
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FIG. 12. Normalized time-headway distribution for the 
TOCA model in free flow and congested traffic for different 
densities. 



The time-headway distributions of the TOCA model, 
however, differ significantly from the results of the NaSch 
model (Fig. 12). For free flow traffic the position of the 
maximum is different from the minimal time-headway for 
the chosen set of parameters. The maximum coincides 
with th while the minimal time-headway is determined 
by the unit of time. For congested traffic the distribution 
has two maxima, one corresponding to the typical time- 
headway in free flow traffic and the other corresponding 
to the typical temporal distance in the outflow region of 
a jam. 

The OV functions of the TOCA and the NaSch model 
differ in two respects (Fig. 13). (i) Due to the fact that 
the randomization step is applied for a finite range of the 
interactions, all cars move deterministically with w max at 
low densities and therefore spatial headways smaller than 
u max cells are completely avoided. This result is at least 
partly a consequence of our simulation setup, i.e., choos- 
ing exactly the same maximal velocity for every car. (ii) 
The second difference is found in the density dependence 
of the OV-function for congested traffic. Because of the 
retarded acceleration in step 1 and the deceleration of 
vehicles with v < d, at very large densities the system 
contains only one large jam with a width comparable to 
the system size. As a consequence, the mean velocity 
at a given distance is reduced considerably compared to 
free flow. The transition to a completely jammed system 
occurs at densities of about 66 veh/km and leads to the 
abrupt change of the OV-curve. 
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FIG. 13. OV functions of the TOCA model in free flow and 
congested traffic for different densities. 



The main difference between the NaSch model and the 
TOCA approach is the structure of jams. Due to the re- 
stricted application of the randomization step, p^ ec must 
be quite large in order to obtain reasonable results for the 
fundamental diagram. Such a choice of Pdco however, re- 
duces significantly the density of jams. This implies that, 
although the typical time-headway in the outflow region 
of a jam has the correct value, the downstream velocity 
of jams is too large. 

Our analysis revealed several shortcomings of the 
TOCA model. However, we believe that the TOCA 
model is an interesting advancement of the NaSch model 
if a finer spatial discretization is applied. We will illus- 
trate this for the example of the density of a wide jam: 
The inverse density of wide jams was used in order to fix 
the size of a cell. This choice is correct, as the jams in 
the model are basically compact, which is not true in case 
of the TOCA model. In this case more accurate results 
could be obtained if each cell would be divided into three 
cells. Using this finer discretization cars occupy two cells 
which would finally lead to a quite realistic dynamics of 
jams. A more elaborate discussion of the discretization 
effects can be found in appendix A. 



IV. CA MODELS WITH MODIFIED DISTANCE 
RULES 

A. The model of Emmerich & Rank 

The CA model introduced by Emmerich and Rank [17] 
(ER-model) is another variant of the NaSch model with 
an enhanced interaction radius. Precisely speaking the 
braking rule of the NaSch model is replaced by applying 
a velocity dependent safety rule that is implemented via 
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a gap- velocity matrix M. The entries of M denote 
the allowed velocities for a car with gap i and velocity 
j. Replacing the braking rule My < j holds because 
otherwise the car would accelerate. For the NaSch model 
the elements of the gap- velocity matrix M^ NaSch ) simply 
read My = min{i, j}. 

Emmerich and Rank tried to improve the NaSch model 
by introducing a larger interaction horizon, i.e. by an 
earlier adaption of the speed. This partly avoids the un- 
realistic effect, that drivers stop from a high speed within 
one time step. Compared to the Nasch model their choice 
of the matrix M only modifies the distance rule for cars 
moving with velocity v max : If 4 < d < 9 the car n has to 
slow down to velocity 4. For all other combinations of d 
and v the NaSch distance rule is left unchanged. 

As a second modification of the NaSch model a differ- 
ent update scheme is applied. The ER model uses an 
unusual variant of the ordered sequential update, i.e., all 
rules, including the movement of the vehicles, are directly 
applied for the chosen car. A unit of time corresponds to 
one update of all cars. Ordered sequential updates use 
normally a fixed sequence of cars or lattice sites. This has 
the disadvantage that some observables, e.g., the typical 
headway, may depend on the position of the detection 
device, even for periodic systems. In order to reduce this 
effect the car with the largest gap is chosen first and than 
the update propagates against the driving direction [17]. 
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FIG. 14. Fundamental diagram of the ER-model. As sug- 
gested in the original work we have chosen 7.5 m as the length 
of a cell, At = 1 s, pace = 0.3 and v max = 5 = 135 km/h. 
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FIG. 15. Normalized time-headway distribution for the ER 
model in free flow and congested traffic for different densities. 

As a consequence of the ordered sequential update 
scheme, the gaps are used very efficiently and very large 
flows can be achieved [43]. (Now, it is allowed that two 
cars arc driving with u max and d — 0, so that flows 
J > 1 veh/At are possible). Therefore large deceleration 
probabilities are necessary to decrease the overall flow to 
realistic values. Nevertheless, due to the sequential up- 
date scheme, the spontaneous jam formation is reduced 
considerably. The application of a sequential update is 
crucial. If it is replaced, e.g., by a parallel update, one 
may observe an unrealistic form, i.e., a non-monotonous 
behavior, of the fundamental diagram at low densities 

[!]• f ^ 

Due to the special choice of M^ ER \ the velocity of 
cars with d < 9 is restricted to v < 4. This means, that a 
generic speed limit with v max = 4 is applied for all den- 
sities p > 1/11 w 12 veh/km, where the mean distance 
between the cars is smaller than 10 cells. Therefore, the 
free flow branch of the fundamental diagram in Fig. 14 
has in contrast to the empirical data two different slopes, 
one corresponding to v max = 5 cells/At if p < 15 veh/km 
and the other to u max = 5 cells/ At at larger densities. 

For the present choice of M we recover basically the 
distance rule of the NaSch model with v max = 4 cells/ At, 
because the speed limit applies only for larger distances. 
Therefore the structure of the congested part of the fun- 
damental diagram is quite similar to the NaSch model. 
However, important differences concerning the micro- 
scopic structure of the traffic state exist, mainly due to 
the modified update scheme. The ordered sequential up- 
date allows motion at high speeds and small distances. 
This could in principle (for small pdoc) lead to very short 
time-headways. For the chosen value of pdec, however, 
the typical time-headways are quite large in the free flow 
regime and do not match the empirical findings. Nev- 
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crthclcss the ordered sequential update changes qualita- 
tively the form of the time- headway distribution, i.e., the 
position of the maximum and the short time cut-off are 
different, as empirically observed (Fig. 15). 

The OV-function of the ER-model differs strongly from 
the empirical findings (Fig. 16). For this quantity the 
modified distance rule is of great importance. In the 
congested regime, we observe plateaus of almost constant 
average velocities v < v max . The density dependence of 
the OV function is, as for the NaSch model, very weak. 
In free flow traffic small headways simply have not been 
observed, in contradiction to the empirical results. 
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FIG. 16. OV function of the ER model in free flow and 
congested traffic for different densities. 



The most important weakness of the ER-model is its 
description of the jam dynamics. First of all for small 
values of pdec the possibility of downstream moving jams 
exist, which contradicts all empirical studies. But even 
for the large value of pdec we applied, jams are not sta- 
ble, i.e., often branch into a number of small jams. There- 
fore it is impossible to reproduce the empirically observed 
parallel moving jams with the ER-model. 

In summary, the gap- velocity matrix allows a more de- 
tailed modeling of the interaction horizon. But keeping 
the parallel update scheme, unrealistic behavior at low 
densities is observed. Using a special variant of the se- 
quential update leads to a very unrealistic structure of 
the microscopic traffic states. 



was introduced in order to provide an alternative mech- 
anism of jam formation. In certain density regimes the 
HS model is very sensitive to external perturbations due 
to its intrinsic nonlinearity. So in contrast to the previ- 
ous approaches the pattern formation is of chaotic rather 
than of stochastic nature, although the definition of the 
model includes a stochastic part as well. 

The deterministic part of the velocity update is done 
by assigning the following velocity to the cars: 



v n {t + 1/2) = v n (t) + X[V op t(d n ) - v n (t)} 



(4) 



where V opt (d) denotes the "optimal" velocity of car n 
for a given distance d n to the vehicle ahead, v n (t) the 
discrete velocity at time t and [• ■ -J the floor function. 
The constant A is a free parameter of the model 5 . The 
acceleration step is the naive discretization of the accel- 
eration step of the space and time continuous OV model. 
In the continuous version of the OV model the parameter 
A determines the timescale of the acceleration. However, 
for time-discrete models it is well known that a simple 
rescaling of time is not possible. Therefore the meaning 
of the parameter A remains unclear. 

The deterministic update is followed by a randomiza- 
tion step as known from the NaSch model, i.e. the veloc- 
ity of a car with v n (t + 1/2) > is reduced with proba- 
bility by one unit. 

Although the definition of the model seems to be quite 
similar to the models discussed in the previous sections, 
important differences exist. In all other models discussed 
so far acceleration is limited to one velocity unit per time- 
step while breaking from w max to zero velocity is possible. 
This is not true for the HS model where a standing car 
may accelerate towards L-^K>pt(oo)J > 1 in a single time- 
step. On the other hand, in particular for small values 
of A, the braking capacity of cars is reduced. A reduced 
braking capacity, however, may lead to accidents (see the 
discussion in [13, 44, 45]), a certainly unwanted feature of 
a traffic model. It also implies that the model is not 
defined completely by the dynamics. This becomes a 
problem especially in simulations. Here further rules are 
necessary to determine how to deal with accidents. 



§In contrast to the original work we consider here only the 
case of one type of cars. Furthermore we denote the OV 
function by V op t so that it can be better distinguished from 
the velocities of the cars. 



B. A discrete optimal velocity model 



Helbing and Schreckenberg (HS) [18] have introduced 
a cellular automata (CA) model for the description of 
highway traffic based on the discretization of the optimal- 
velocity (OV) model of Bando et al. [33]. The model 
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FIG. 17. Fundamental diagram of the HS-model. As sug- 
gested in the original work we have chosen 2.5 m as the 
length of a cell, a vehicle has a length of 2 cells, At = 1 
s, Pdcc = 0.001, Umax = 15 = 135 km/h and A = 1/1.3. 
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FIG. 18. Normalized time-headway distribution for the HS 
model in free flow and congested traffic for different densities. 



We will discuss the possibility of accidents in some 
more detail in appendix B. This discussion concentrates 
on a criterion which ensures that for any possible initial 
condition no accident occurs. 

In [18] for comparison with empirical data the following 
OV-function is suggested: 
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The length of a cell is set to Ax = 2.5 m, At = 1 is 
chosen as the unit of time, A = j-^ and we used the ran- 
domization probability pdoc = 0.001 as suggested in [18]. 
A vehicle has a length of I = 2 cells corresponding to 
5 m. For this choice of A and the OV function the model 
is not strictly free of collisions as our discussion in the 
appendix shows, but does at the same time not lead to 
accidents if an appropriate initial condition is chosen and 
the density is not too high. 
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FIG. 19. Distance headway distributions in the congested 
regime. Obviously non-compact jams coexist with free flow 
regimes, where the distance between two cars is rather large. 
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FIG. 20. OV function of the HS model in free flow and 
congested traffic for different densities. 

The optimal velocity function that governs the deter- 
ministic part of the vehicle dynamics, leads to speed 
limits in certain density regimes Z/(d max + V) < p < 
l/(dmin+l) with V r op t(rfmin) = V opt (d max ). These different 
optimal speeds become visible in different slopes in the 
free flow branch of the fundamental diagram (Fig. 17). 
For congested traffic two different traffic regimes can be 
identified, as empirically observed. For very high densi- 
ties, one observes a reasonable agreement with the em- 
pirical data, i.e. the form of the jammed branch is re- 
produced qualitatively. This branch of the fundamental 
diagram is, however, observed only in a very narrow in- 
terval of global densities. 

Compared to the two other traffic states the repro- 
duction of synchronized traffic is rather poor. First, one 
obviously observes a strong correlation between density 
and flow, which is contrast to the empirical findings, and 
second, the range of densities which is observed in local 
measurements is quite narrow. 

The main difficulties of the model are visible when 
comparing it with empirical results on a microscopic 
level. The simulations for the time-headway distribu- 
tion show a strong density dependence of the maximum 
for the free flow states. This is due to the long-ranged 
interactions that tend to generate traffic states that are 
very homogeneous. Therefore short time-headways are 
suppressed at low densities. The second problem is the 
quasi-deterministic character of the model. This implies 
that drivers obey the distance rule in almost any case. 
As a result the peak values of the time-headway distribu- 
tion have extremely high weights. In congested traffic wc 
observe a density independent position of the maximum 
of the time-headway distribution. The maximum carries 
almost the whole weight of the distribution, in contradic- 



tion to the empirical findings. The reason for this can be 
read off from the distance headway distributions for dif- 
ferent global densities (Fig. 19). Within a large density 
regime we observe coexistence of non-compact jams and 
free flow traffic. Therefore we can state that both high 
density states correspond to stop-and-go traffic, i.e. the 
model fails to reproduce synchronised traffic at all. 

The mismatch of the model and empirical structure of 
traffic states is also obvious for the OV-function (Fig. 20). 
It shows almost no density dependence and is basically 
independent of the traffic state. The difference between 
the different curves is only in a density dependent cut-off 
of the distribution, i.e. at high densities large distances 
simply do not occur. 

The simulations show that HS model fails to reproduce 
the microscopic structure of the empirical observed traffic 
states. From our point of view the problems describing 
the empirical observation are due to the nature of the 
model. It introduces a static rule that leads to a reason- 
able agreement with the empirical fundamental diagram. 
For a proper choice of A the vehicles take instantaneously 
a velocity close to the optimal velocity, i.e. the dynam- 
ical aspects of highway traffic are extremely simplified. 
Therefore inhomogencous traffic states are only observed 
in the presence of quenched disorder [18], e.g. different 
types of cars, and not produced spontaneously. 

V. BRAKE LIGHT VERSION OF THE NASCH 
MODEL 

Quite recently a brake light (BL) version of the NaSch 
model has been introduced [19, 20] in order to give a 
more complete description of the empirically observed 
phenomena in highway traffic. In contrast to the models 
we considered in the previous sections, which represent 
already well known modeling approaches, we also discuss 
the basic features of the model that have not been pre- 
sented so far. In the development of the model the main 
aim was the reproduction of the empirical microscopic 
data in a robust way. 

A. Definition of the BL model 

The BL model combines several elements of older mod- 
eling approaches, e.g., velocity anticipation [46,47] and a 
slow-to-start rule [15,39]. In addition, a dynamical long 
ranged interaction is included: In their velocity depen- 
dent interaction horizon drivers react on brakings of the 
leading vehicle that arc indicated by an activated brake 
light [48]. The interaction, however, is limited to nearest 
neighbor vehicles [49]. The update rules are formulated 
in analogy to the VDR model. In particular the inter- 
actions are strictly local and a parallel update scheme is 
applied. 
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In order to allow for a finer spatial discretization for a 
given length of a car, we include the possibility that a car 
may occupy more than a single cell. Therefore the gap 
between consecutive cars is given by d n = x n+ \ — x n — I 
(where I is the length of the cars). The brake light b n can 
take on two states, i.e., on (off) indicated by b n = 1(0). 
In our approach the randomization parameter pdec for 
the nth car can take on three different values po, Pd and 
Pb, depending on its current velocity v n (t) and the status 
b n+ \ of the brake light of the preceding vehicle n + 1: 

Pdcc = Pdee(v n (t),b n+ i(i),th,t s ) 




p b if b n+1 = 1 and t h < t, 

Po if v n = 

Pd in all other cases. 



The two times th = -^jt) anc ^ *s = m ^ n l v n(t) , h} , 
where h determines the range of interaction with the 
brake light, are the time th needed to reach the position 
of the leading vehicle which has to be compared with 
a velocity-dependent (temporal) interaction horizon t s . 
t s introduces a cutoff that prevents drivers from react- 
ing to the brake light of a predecessor which is very far 
away. Finally d!° ff) = d n + max{w anti - <4ecurity, 0} de- 
notes the effective gap where v an ti = min{d„+i, is 
the expected velocity of the leading vehicle in the next 
time-step. The effectiveness of the anticipation is con- 
trolled by the parameter ^security Accidents are avoided 
only if the constraint ^security > 1 is fulfilled. The update 
rules then are as follows (t < t\ < t 2 <t+l): 

0. Determination of the randomization parameter: 

Pdcc = Pd C c(v n (t),b n+ i(t),t h ,t s ) 

b n {t + l) = 

1. Acceleration: 

if ((6n+i(t) = 0) and (b n (t) = 0)) or (t h > t s ) then 
v n (h) = mm{v n (t) + l,v max } 

2. Braking rule: 

v n (t 2 ) = min{dn S \v„(ti)} 
if (v n (h) < v n (t)) then 

M* + i) = i 



Here rand() denotes a uniformly distributed random 
number from the interval [0, 1]. 

The new velocity of the vehicles is determined by steps 
1 — 3, while step determines the dynamical parameters 
of the model. Finally, the position of the car is shifted in 
accordance with the calculated velocity in step 4. 

In order to illustrate the details of the approach we 
now discuss the update rules step- wise. 

0) The braking parameter pdcc is calculated. For a 
stopped car the value pace = Po is applied. There- 
fore po determines the upstream velocity of the 
downstream front of a jam. 

If the brake light of the car in front is switched 
on and it is found within the interaction horizon 
Pdcc — Pb is chosen. A car perceives a brake light 
of the vehicle ahead within a time dependent in- 
teraction horizon t s = min{w„(i), h} where v n (t) 
is the current velocity and h an integer constant. 
The velocity dependence takes into account the in- 
creased attention of the driver at large and reduces 
the braking readiness at small velocities. This reac- 
tion is performed only with a certain probability of 
Pb- In order to obtain a finite range of interactions 
a cutoff at a horizon of h seconds is made **. 

Finally, pdcc = Pd is chosen in all other cases. 

1) The velocity of the car is increased by one unit (if 
it does not already move with maximum velocity). 
The car does not accelerate if its own brake light or 
that of its predecessor is on and the next car ahead 
is within the interaction horizon. 

2) The velocity of the car is adjusted according to the 
effective gap. 

The brake light of a vehicle is activated only if the 
new velocity is reduced compared to the preceding 
time-step. Note that the application of the braking 
rule does not necessarily lead to a change of the ve- 
locity, as it can compensate a previous acceleration. 
The restriction stabilizes dense traffic flows. 

3) The velocity of the car is reduced by one unit with 
a certain probability p dcc = Pd cc {v n {t), 

If the car brakes due to the predecessor's brake 



3. Randomization, brake: 
if (rand() < pa cc ) then { 

v n (t + 1) = iasx{v n (t 2 ) - 1, 0} 
if ((Pdcc = Pb) and (v n (t + l) = v n (t 2 ) - 1)) then 
b n (t + l) = 1} 

4. Car motion: 

(t + l) = x n (t) +v n (t + 1) 



"Indeed, increasing pdcc to pb is the simplest possible re- 
sponse to the stimulus brake light. More sophisticated re- 
sponse functions like a direct reduction of the velocity or the 
gap are conceivable but lead to some problems in combina- 
tion with anticipation. In addition, one can think of different 
implementations of the brake noise p b . For example, we have 
tried more sophisticated p(,-functions, like a linear relation- 
ship between p b and the velocity, the difference velocity to 
the predecessor or the gap, but for the sake of simplicity in 
this paper we will focus on a constant p b . 
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light, its own brake light is switched on. We also 
stress the fact that even for distances d n < h the 
action of the brake light is restricted to brakings 
that are induced by the vehicle in front (either by 
the braking rule or by an activated brake light) and 
not by spontaneous velocity fluctuations. 

4) The position of the car is updated. 



car in front is indeed a strong stimulus to adjust the own 
speed. Therefore pb has typically a high value. Finally, 
^security tunes the degree of the velocity anticipation and 
has a strong influence on the cut-off of the time-headway 
distribution. 



C. Validation of the full model 



B. Calibration of the model 

The following parameters of the model allow to ad- 
just the simulation data to the empirical findings: the 
maximal velocity w max , the car length I, the braking pa- 
rameters pd, Pb, Po, the cut-off h of interactions, and the 
minimal security gap rf sccur i ty . The parameters of the 
have been chosen such that they can easily be related to 
the empirical findings. As in the previous models a single 
set of parameters is used for all traffic states. 

In order to obtain realistic values of the acceleration 
behavior of a vehicle, the cell length of the standard CA 
model is reduced to a length I of 1.5 m. Since the time- 
step is kept fixed at a value of 1 s this leads to a veloc- 
ity discretization of 1.5 m/s which is of the same order 
as the "comfortable" acceleration of somewhere about 
1 m/s 2 [50]. Like in the standard CA model a vehicle 
has a length of 7.5 m that corresponds to 5 cells at the 
given discretization (see appendix A for a discussion of 
the discretization effects). 

Some of the parameters can be fixed as, e.g., in the 
VDR-model: The maximum velocity w m ax is determined 
by the slope of the free flow branch of the fundamental 
diagram. The upstream velocity of a jam can be tuned by 
the parameter po and the strength of fluctuations that are 
controlled by the parameter p d determine the maximal 
flow. 

The other parameters of the model are connected with 
an interaction that have not been included in the models 
we discussed so far. The parameter h describes the hori- 
zon above which driving is not influenced by the leading 
vehicle. Several empirical studies reveal that h corre- 
sponds to a temporal headway rather than to a spatial 
one. The estimates for h vary from 6 s [51], 8 s [52,53], 
9 s [54] to 11 s [55]. Another estimation for h can be ob- 
tained from the analysis of the perception sight distance. 
The perception sight distance is based on the first per- 
ception of an object in the visual field at which the driver 
perceives movement (angular velocity). In [56] velocity- 
dependent perception sight distances are presented that, 
for velocities up to 128 km/h, are larger than 9 s. We 
therefore have chosen h to be 6 s as a lower bound for 
the time-headway. Besides, our simulations show, that a 
good agreement with empirical data can only be obtained 
for h > 6. This corresponds to a maximum horizon of 
6 x 20 cells or a distance of 180 m at velocity v max . 

The next parameter one has to fix is pb- This parame- 
ter controls the propagation of the brake light. A braking 



With this parameter set we have calibrated the model 
to the empirical data. Leaving p , h and v max fixed, 
we got the best agreement with the empirical data for 

Pdoc = 0.1, Pb = 0.94 and Security = 7. 

As one can see in Fig. 21 the slope of the free flow 
branch and the maximum flow coincides with the empir- 
ical data indicating that v max and pdec have been chosen 
properly 
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FIG. 21. Local fundamental diagram obtained by the sim- 
ulation of the brake light version of the NaSch model. The 
parameters are: po = 0.5, h = 6, « m ax = 20, Pdec = 0.1, 
Pb — 0.94, d sccur ity = 7. A time-step corresponds to At — 1 s, 
a cell has a length I — 1.5 m and a vehicle covers 5 cells. 



However, the simulated densities are less distributed 
than in the empirical data set. The width of the density 
distribution is of the same order as it was found for the 
NaSch and ER-model. The mismatch between simulation 
and empirical results of the density can be related to 
discretization effects, which introduce an upper limit for 
the density if simple (virtual) counting loops are used as 
detection devices. 

A second lower branch appears for small values of the 
flow which represents wide jams. Because only moving 
cars are measured by the inductive loop large densities 
cannot be calculated, as in the empirical data of Fig. 1. 
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FIG. 22. Cross-correlation of the flow and the density in 
free flow and congested traffic for different densities and ho- 
mogeneous initialization. 

The next parameter that can be directly related to 
an empirical observable quantity, namely the upstream 
velocity of the downstream front of a wide jam, is the 
deceleration probability pq. 

We used the calculation of the density autocorrelation 
function in the congested state of a system that was ini- 
tialized with a mega jam for the determination of the 
velocity of the jam front. One obtains an average jam 
velocity of 2.36 cells/s (=12.75 km/h) for p a = 0.5. This 
jam velocity is independent of the traffic condition and 
holds for all densities in the congested regime. Thus, al- 
though metastable traffic states can be achieved by the 
finer discretization (see appendix A) the slow-to-start 
rule is necessary for the reduction of the jam velocity 
from about 20.45 km/h to 12.75 km/h. This velocity is 
also in accordance with empirical results [32] . 

In Fig. 22 the cross-covariance cc( J, p) of the flow and 
the local measured density for different traffic states is 
shown. In the free flow regime the flow is strongly cou- 
pled to the density indicating that the average velocity is 
nearly constant. Also for large densities, when wide jams 
are measured, the flow is mainly controlled by density 
fluctuations. In the mean density region there is a tran- 
sition between these two regimes. At cross-covariances 
in the vicinity of zero the fundamental diagram shows a 
plateau. Traffic states with cc(J, p) w were identified 
as synchronized flow [22]. In the further comparison of 
our simulation with the corresponding empirical data we 
used these traffic states for synchronized flow data and 
congested states with cc(J, p) > 0.7 for data of wide jams. 
The results show that the approach leads to realistic re- 
sults for the fundamental diagram and that the model is 
able to reproduce the three different traffic states. 

To characterize the three traffic states, we calculated 



the autocorrelation of the flow, the density as well as the 
velocity for different global densities. 

In free flow, the density and the flow show the same 
oscillations of the autocorrelation function, whereas the 
speed is not correlated in time. 

In contrast to the NaSch model, the autocorrelation 
function at large densities shows a strong coupling of the 
flow and the velocity. Now, the velocity of a car not only 
depends on the gap but also on the density, so that the 
flow and the velocity are mainly controlled by the density 
(Fig. 23). 




T [min] 

FIG. 23. Autocorrelation function of the density, the ve- 
locity and the flow for p — 67 veh/km with a random initial- 
ization. 

Next we compare the empirical data and simulation 
results on a microscopic level. 

In Fig. 24 the simulated time-headway distributions 
for different density regimes are shown. 

Due to the discrete nature of the model, large fluctua- 
tions occur and the continuous part of the empirical dis- 
tribution shows a peaked structure at integer-numbered 
headways for the simulations. In the free flow state ex- 
tremely small time-headways have been found, in accor- 
dance with the empirical results. This is qualitatively 
different from the other CA models with parallel update 
scheme. 
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boundary conditions, we observe a smooth time head- 
way distribution, which is in good agreement with the 
empirical data (see Fig. 25). Therefore we can state that 
properties like the width as well as the smoothness of 
the time-headway distribution are strongly dependent on 
the choice of the simulation setup. In contrast, the short- 
time cut-off of the distribution is model dependent. Time 
headways shorter than the chosen unit of time are in case 
of a parallel update only observed if anticipation effects 
are included. The actual value of the cut-off for a given 
unit of time is tuned by the parameter d sccur it y . The re- 
sults for congested flow, however, are not influenced by 
different types of vehicles. 

The ability to anticipate the predecessor's behavior be- 
comes weaker with increasing density so that the weight 
of the small time- headways is reduced considerably in the 
synchronized state. The maximum of the distribution 
can be found in the vicinity of 1 s in accordance with the 
empirical data, the density dependence, however, cannot 



FIG. 24. Time-headway distribution for different densities 
in free flow (top) and in the synchronized state (bottom) . 




FIG. 25. Time-headway distribution in the free flow regime 
of a system with open boundary conditions and different types 
of vehicles. The maximal velocity of the slow vehicles was set 
to as «max = 108 km/h = 20 cells/s and of the fast vehicles 
as «max = 135 km/h = 25 cells/s. We considered 15% of 
the vehicles as fast vehicles (note that these are vehicles that 
disregard the speed limit). 



Nevertheless, for our standard simulation setup at 
small densities the statistical weight of these small time- 
headways is significantly underestimated. This apparent 
failure of the model is the result of the chosen simulation 
setup. If we introduce different types of cars and open 
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FIG. 26. The OV-function for different densities in the free 
flow and congested regime. 

Instead, with increasing density the maximum at a 
time of 1 s (in the NaSch model the minimal time- 
headway is restricted to 1 s because of rule 2) becomes 
more pronounced. This result is also due to the dis- 
cretization of the model that triggers the spatial and 
temporal distance between the cars. Because of the ex- 
ponential decay of the waiting time distribution of cars 
leaving a jam, the peak at a time of 1 s is the most prob- 
able in the time-headway distribution. 

The OV curve of our model approach shows an ex- 
cellent agreement with empirical findings. For densities 
in the free flow regime it is obvious that the OV-curve 
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(Fig. 26) deviates from the linear velocity-headway curve 
of the NaSch model. Due to anticipation effects, smaller 
distances occur, so that driving with w max is possible even 
within very small headways. This strong anticipation be- 
comes weaker with increasing density and cars tend to 
have smaller velocities than the headway allows so that 
the OV-curve saturates for large distances. 
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FIG. 27. Autocorrelation of the speed and of the spatial 
and temporal headway for free flowing vehicles (top) and 
for a synchronized state (bottom). In order to obtain a 
slow decay of the speed autocorrelation function in the free 
flow regime the simulation was performed on a open sys- 
tem with 20% of slow cars (v s ^°™ = 15 cells/s = 81 km/h, 
Vmfx = 20 cells/s = 108 km/h.) 

The saturation of the velocity, which is characteristic 
for synchronized traffic, was not observed in earlier ap- 
proaches. The value of the asymptotic velocities can be 
adjusted by the last free parameter pf,. The OV-curve 
in the synchronized regime is independent of the maxi- 
mum velocity and is only determined by the dynamical 
behavior of the model. 

Next, we calculated the autocorrelation of the time- 
series of the single- vehicle data (Fig. 27). Note, that the 
data of the free flow state was collected in an open sys- 
tem with 20% of slow cars with v max — 15 cells/s = 81 
km/h. In the free flow regime the data shows a strong 
coupling of the spatial and temporal headway that sup- 
ports the results obtained by aggregated data (J oc At -1 
and p oc Ax -1 ). In contrast, the autocorrelation of the 
velocity shows a slow asymptotic decay. This supports 
the explanation of [22] that the slow decrease for small 
distances is due to small platoons of fast cars led by one 
slow car. In the synchronized state, longer correlations 
of the speed and the spatial headways can be observed. 
So, similar to the free flow regime, in the synchronized 
regime platoons of cars appear that are moving with the 



same speed [24]. 

VI. THE MODEL OF KERNER, KLENOV AND 
WOLF 

The most recent modeling approach we include in our 
comparison was introduced by Kerner, Klcnov and Wolf 
(KKW) [21]. This model is a fully discretized version 
of the space-continuous microscopic model introduced 
by Kerner and Klenov [58]. It combines, as the BL 
model, elements of car-following theory with the stan- 
dard distance-dependent interactions. It is defined by 
an update rule including a deterministic and a stochastic 
part. The deterministic rule, which reads (t < t\ <t+l) 

u„(*i) = max|o,min{w max , v sa f c (t),-y dos (t)}|, (6) 

is applied first. The three velocities appearing are the free 
flow or maximal speed of the cars v max , the safe velocity 
^safe(t) an d finally the desired velocity Vd es (t). v sa f c (t) is 
the velocity which guarantees collision-free motion and is 
simply the gap to the preceeding car, v aa ,i c (t) = d n (t). It 
is the introduction of t>des(i) which makes the difference 
to the NaSch model. The velocity Vd es (t) is given by 

, . _ J v n (t) + a for d n > D(v n (t)) - I, . . 

Vdes[t) ~ \ v n (t) + A(t) for d n < D(v n (t)) - I. (l) 

The calculation of Vdes(t) replaces the acceleration step 
of the NaSch model by a more complex rule. Here I is 
the length of the vehicles and D(v) a synchronization 
distance. The authors suggested a linear 

D(v) = D a + kv. (8) 

and a quadratic form 

D(v) =D Q + v + (3v 2 (9) 

for the velocity dependent interaction range. Apart from 
the task of choosing an appropriate function, two model 
parameters are introduced in both cases. So far no sys- 
tematic analysis of traffic data exist which leads empir- 
ically based parameter values or the functional forms. 
This could be done, at least in principle, by an exten- 
sive analysis of floating-car measurements. Moreover, 
the results in [21] show that the results agree at least 
qualitatively for both functions (8), (9) which have been 
considered. 

The interaction range has been introduced as a syn- 
chronization radius, i.e. D(y) is the distance which sep- 
arates free driving cars from cars which already adjust 
their velocity according to the vehicle ahead. For large 
distances to the vehicle ahead, d n > D(v n (t)) — l, the cal- 
culation of fdos is equivalent to the acceleration step of 
the NaSch model. Inside the enlarged interaction radius, 
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however, v^cs depends on the velocity of the leading car. 
Explicitely A(t) is given by 

{-b if v n (t) > v n+ i(t) 
if v n (t) = v n+1 (t) (10) 
a if v n (t) < v n+1 (t), 

This means that within the interaction radius drivers 
tend to adapt their velocity to the vehicle ahead. 
The second update rule is stochastic. It is given by 



The stochasticity is included in the term v n (ti) + 
while the others are in order to guarantee that the new 
velocity is below the speed limit, leads to no collisions and 
is in accordance with the chosen acceleration capacity 
a n of the cars. The stochastic variable r] can take the 
following values: 

{-1 if r < p b , 
1 if p h < r < pb+p a , (12) 
otherwise 

Both probabilities p a and pb introduced here are velocity 
dependent. One has 

, x f po if v = , 10 s 
P^ = \p i{v> 0. ( 13 ) 

with po > p. The stochastic braking is analogous to the 
slow-to-start rule known from the VDR model. Contrary 
the stochastic acceleration is a new feature of the model 
which weakens the synchronization of speeds as it applies 
to cars which reduced or kept their velocity although safe 
driving would have allowed a larger velocity. The func- 
tion p a (vn) is explicitly given by 

PM = \ P&1 l V V P ( 14 ) 

V ' \ Pa2 if V > V p , V ; 

where v p , p al and p a2 < p a \ arc adjustable parameters of 
the model. The different probabilities have to be chosen 
such that p a + Pb < 1 is fulfilled for any velocity. The 
velocity update is completed by this second stochastic 
rule and is followed by a parallel update of the positions. 

For further illustration of the update rules we compare 
them briefly to the BL model. Both models include the 
update rules of the VDR model and enlarge the inter- 
action radius of the drivers within a velocity dependent 
interaction range. The driving strategy within this larger 
interaction range is, however, different. While the BL 
model introduces an event driven interaction model, the 
KKW is more car-following like. Another important dif- 
ference is that the velocity anticipation is not included 
in the approach of [21], although such an extension is 
possible [59]. 
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FIG. 28. Local fundamental diagram of the KKW model 
for the following set of parameters: The length of a cell is 
set to 0.5 m. Each car occupies I — 15 cells. The maximal 
velocity is given by = 108 km/h = 60 cells/At, where 

At = 1 s. Also the other model parameters are set to the 
values suggested in [21]: a = b = 1, Do = 60, k = 2.55. The 
parameters determining the stochastic part of the model take 
the values: p = 0.04, p = 0.425, p a i = 0.2, p a2 = 0.052 and 
v p = 28. 

Fig. 28 shows the fundamental diagram of the KKW 
model, obtained by local measurements flux and density 
in a periodic system. Compared to the other models 
we analyzed one observes two remarkable differences: In 
synchronized traffic the flow has a local minimum for a 
density of 30 veh/km and reaches a second maximum 
for a density of 40 veh/km. The origin of this structure 
lies in the stochastic acceleration of cars which reduces 
considerably the probability to form a jam. A second 
important feature is the complex structure of the fun- 
damental in the presence of jams. For very high global 
densities one observes all three traffic states at the same 
time and no strict phase separation as, e.g., in case of the 
VDR model. 



v n (t + l) = max{0,min{w„(ii) + (tl) + a, Vf roo , »max}}, 

(11) 
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FIG. 29. OV function in free flow and congested traffic 
of the KKW model for different densities. The same set of 
parameters as in Fig. 28 has been used. 



Measurements of the OV-function show that the mi- 
croscopic structure of the model differs from the empir- 
ical findings. In free flow traffic small headways are al- 
most not observed and the maximum speed is reached at 
larger distances than in real traffic. This indicates that, 
compared to real traffic, the repulsive part of the car-car 
interactions is overemphasized. While the differences be- 
tween empirical data and model results might be reduced 
for a different set of model parameters, the model results 
for synchronized traffic differ even qualitatively. In real 
traffic one observes for a given density a crossover from 
a density independent form of the OV-function at small 
distances to an asymptotic velocity for larger distances 
(see Fig. 4) which depends on the density. This is not 
reproduced by the KKW model, where a distance inde- 
pendent average velocity is observed only in a narrow 
range of spatial headways, if it is observed at all. 

The comparison between empirical and simulation re- 
sults of the time-headway distribution indicates that the 
model largely fails to reproduce the empirical results ob- 
tained for free flow traffic. This is, as discussed before, 
partly a results of the simplified setup we used for our 
simulations. A much better agreement would be obtained 
if we consider a realistic distribution of maximal veloc- 
ities. But even in this case one is left with a problem. 
The lack of velocity-anticipation leads to a sharp cut- 
off of the time-headway distributions for times less than 
one unit of time, i.e. 1 s. Although the position of the 
cut-off can be tuned by varying the temporal discretiza- 
tion, it must be noted that this still does not lead to 
the right functional form, as the maximum of the time- 
headway distribution is located at the minimal observed 
time-headway. This again confirms necessity of velocity 



anticipation for the reproduction of the empirical findings 
at short time-headways. 

Summarizing the CA model introduced by Kerner, 
Klenov and Wolf reveals three distinguishable traffic 
states, as observed in empirical studies. The reproduc- 
tion of the empirical time-headway distribution and fun- 
damental diagram is partly satisfying and could be eas- 
ily improved by the introduction of velocity anticipation. 
The most important differences between empirical find- 
ings and model results concern the OV-function. This in- 
dicates that the microscopic structure of the model states 
does not match the real structure of highway traffic. We 
also believe that this disagreement is due to the very na- 
ture of the car-car interactions in the KKW model and 
cannot be resolved by a better choice of the model pa- 
rameters. 



VII. COMPARISON OF THE FUNDAMENTAL 
DIAGRAMS 

The comparison of the models presented so far is based 
on local measurements of inductive loops. Therefore, the 
model parameters have been chosen in order to allow the 
best possible accordance with the empirical setup. One 
of the main disadvantages of local measurements is that 
the detected values of the flow and the velocity strongly 
fluctuate whereas density cannot even be defined locally 
in a strict sense. However, in traffic flow simulations it is 
possible to get averaged quantities that are representa- 
tive for a given density. Therefore, in this section global 
measurements of the flow and the density of the various 
models are given for a typical set of parameters in or- 
der to demonstrate the characteristics of the approaches. 
However, since density can be calculated exactly, the dis- 
tinction between the traffic states is omitted. 
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FIG. 30. Comparison of the global fundamental diagram 
of the NaSch model, the VDR model, the TOCA, and the 
ER model for typical parameter values and a homogeneous 
initialization. 



Density, flow and velocity can be measured globally 
in the following way: The density p g i bai can directly 
be obtained by counting the number TV of vehicles on a 
highway section of length L via 

Pglobal = (15) 

The average velocity v g i b a i is then defined as 

1 N 

^global = Vn ( 16 ) 

n— 1 

with the velocity v n of vehicle n. Again, the hydrody- 
namical relation allows the calculation of the flow 

1 N 

^/global = PglobalWglobal = ^ V n . (17) 

n— 1 

A typical fundamental diagram consists of a linear free 
flow branch that intersects with an almost linear con- 
gested branch. As one can see in Fig. 30 and Fig. 31 
nearly all discussed models are able to reproduce this 
basic characteristics. The fundamental diagram of the 
HS model, however, exhibits two distinct maxima. The 
first maximum is simply given by the transition from free 
flow to congested traffic. The second maximum is a con- 
sequence of the chosen OV-curve. Since vehicles with 
3 < d < 5 have to drive with a velocity of 2, the flow in- 
creases linearly for densities well above a certain density 
until the average gap is smaller than 3 cells. Moreover, 
due to the OV-curve the vehicles behave deterministi- 
cally and choose their velocity according to the gap, i.e. 
v = d. As a result of a nearly uniform gap distribution, 
effectively speed limits are applied for certain density in- 
tervals which are reflected by the occurrence of different 
slopes in the congested branch of the fundamental dia- 
gram. 
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FIG. 31. Comparison of the global fundamental diagrams 
of the HS model, the BL model, and the KKW model with 
that of the NaSch model for typical parameter values and a 
homogeneous initialization. 

This behavior is typical for models with modified dis- 
tance rules and can also be found in the ER model. Since 
the choice of the gap-velocity matrix in the ER model 
leads to speed limits for different density regimes, the 
free flow branch shows two different slopes like in the 
local measurements. Even more severe is the lack of a 
distinct maximum in the fundamental diagram. This is 
a consequence of the ordered sequential update of the ER 
model. It is possible that jams can also move in down- 
stream direction, thus leading to many small jams with 
a large flow. 

Measurements of empirical data have revealed that the 
outflow from a jam is reduced considerably compared to 
the maximum possible flow. As a result, metastable free 
flow states exist and hysteresis effects can be observed in 
the fundamental diagram [57]. 

Obviously, this is the case for the VDR model, the 
TOCA model as well as for the BL and KKW models 
while the maximum possible flow of the NaSch model is 
as large as the outflow from a jam. 

Since the deceleration probability in the VDR model 
was chosen very small (pdec = 0.01) the stability of the 
homogeneous branch of the fundamental diagram is very 
large. In contrast, once a jam has formed above a certain 
threshold density the large deceleration probability for 
the vehicles at rest is responsible for the reduced outflow 
from a jam. As a result, the system is phase separated 
into a region of free flow and a compact moving jam. 
The capacity drop can simply be tuned by varying the 
difference between the two deceleration parameters. In 
analogy to the VDR model, in the TOCA model only ve- 
hicles with v < d decelerate with the probability p. Thus, 



22 



vehicles driving with v max and d > v lead to a stable high 
flow branch in the fundamental diagram up to a density 
of - — ~fi- However, the congested regime of the TOCA 
model reveals the existence of two different slopes in the 
fundamental diagram. For densities larger than 1/2 ve- 
hicles have on average a gap of less than one cell. Since 
the vehicles decelerate with a large probability, but do 
accelerate with a rate smaller than one, the system now 
contains only one large jam whose width is comparable 
to the system size. 
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FIG. 32. Fundamental diagram for different horizons h 
(top) and for different pb (bottom). All simulations have been 
performed with an homogeneous initialization. 



Like in the VDR model, in the BL model the high 
flow states can simply be controlled by the deceleration 
parameter p for vehicles at rest. However, in the con- 
gested regime two distinct slopes of the fundamental di- 
agram become visible. The density at which the slope 
changes and the shape of the fundamental diagram can 
be triggered by the parameters h and pb that determine 
the interaction between vehicles with d > v. In partic- 
ular, the higher h the smaller the density p max of the 
maximum flow (Fig. 32 top). For large h the fundamen- 
tal diagram converges very fast, so that the fundamental 
diagram for values larger than h = 8 are identical. More- 
over, even small values of pb have a strong influence on 
the flow. The high flow branch of the fundamental dia- 
gram (Fig. 32 bottom) and the density p max of maximum 
flow are reduced. For large values the congested branch 
of the fundamental diagram shows two different slopes. 
The higher pb, the smaller the density at which the slope 
changes. 



A. Minimal model? 

The BL model improves, compared to the other mod- 
els we discussed in this work, the agreement with the 
empirical data, especially in the case of the OV-curve. 
Nevertheless, this is only possible with the application 
of a variety of new update rules. Therefore, it remains 
an open question whether this set of update rules can be 
reduced. 

In the top part of Fig. 33 we successively dropped the 
extensions of the model. First, the slow-to-start rule has 
been omitted. Without the slow-to-start rule the model 
lacks the ability of a reduced outflow from a jam and 
the number of large compact jams is reduced so that the 
flow increases. As a further reduction of the model, antic- 
ipation is switched off. This leads to a decrement of the 
flow at densities larger than the density of maximum flow. 
Now headways smaller than the velocity are not possible, 
which manifests in the OV-curve at small densities. For 
large densities the anticipation of the predecessors veloc- 
ity becomes more and more difficult until anticipation is 
no longer applicable. Therefore, the differences between 
the curves with and without anticipation vanishes. 
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FIG. 33. Successive extension of the NaSch model with 
brake lights (top) and without brake lights (bottom). Note 
that the system is initialized in an homogeneous state to gen- 
erate also high flow states. 



Applying the braking rule as the only extension leads 
to a plateau-like fundamental diagram compared to the 
NaSch model. Additionally, the flow is reduced dramat- 
ically. It is the brake rule that changes the shape of the 
fundamental diagram. 
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In bottom part of Fig. 33 the same successive reduc- 
tions of the rules have been applied to the model without 
braking rule. Neither the anticipation, nor the slow-to- 
start rule applied as a single extension or in combination 
are able to change the shape of the fundamental diagram. 

Considering the empirical fact that small time- 
headways and a reduced outflow from a jam exist, the 
braking rule is the only new extension of the NaSch 
model. This new rule turns out to be crucial for the 
correct generation of the OV-curves and the occurrence 
of synchronized traffic. 

So the set of rules chosen for the BL model is minimal 
in the sense that all are needed to obtain a satisfactory 
agreement with empirical data. We also believe that is 
essential to combine car-following like behavior and dis- 
tance based rules. In case of the BL-model the velocity 
adjustment is event driven, i.e. the drivers react to brak- 
ing cars in the upstream flow. It is not excluded that the 
same can be achieved with a different, but simpler set 
of rules. This is highly desirable in order to reduce the 
complexity of the model and the number of parameters. 
However, it is currently unclear whether there is a simi- 
larly simple physical mechanism behind the formation of 
synchronized traffic as it is behind the formation of wide 
jams. For the latter, the reduction of the outflow from 
a jam below the maximal flow is essential which can be 
easily achieved by any kind of slow-to-start rule. 

VIII. DISCUSSION 

The intention of our investigation was to single out 
the models which are able to describe the empirically ob- 
served microscopic structure of traffic flow correctly. It is 
well-known that many quite different models exist which 
reproduce the macroscopic properties (e.g. global funda- 
mental diagrams or spontaneous jam formation) rather 
accurately [1-3]. However, recently single- vehicle data 
have become available. A thorough analysis of these data 
has allowed for a deeper understanding of the microscopic 
properties which now should be incorporated into the dif- 
ferent modeling approaches. 

We have suggested a test scenario based on the com- 
parison of computer simulations in a realistic setup with 
empirical data obtained using stationary inductive loops. 
An important point is that we have used only one fixed 
set of model parameters which has been determined by 
comparison with empirical data, e.g. with the free-flow 
velocity. Therefore we are able to determine whether a 
model is able to describe all traffic situations consistently 
without the necessity to tunc parameters according to the 
state. 

Our focus was on cellular automata models [1] and es- 
pecially variants of the Nagel-Schreckenberg model [10, 
11] which can be considered as a minimal CA model for 
traffic flow. Our comparison has revealed differences be- 
tween the models on a macroscopic scale which become 



even more pronounced on a microscopic level of descrip- 
tion. 

We have seen that models with modified distance rules, 
like the ER and the HS model, have problems on a macro- 
scopic level. They are not able to produce a realistic 
(global) fundamental diagram and it is difficult to make 
these models intrinsically crash-free. 

The NaSch model, the VDR model, the TOCA model 
and the brake light version of the NaSch model reproduce 
the fundamental diagram quite well. This is already suffi- 
cient for many applications. In urban traffic, for example, 
the dynamics of the vehicles between two intersections is 
predominantly determined by traffic lights. The correct 
description of queues at cross-roads, therefore, only re- 
quires the existence of two distinct traffic phases, namely 
free flow and congested traffic. 

More realistic applications of traffic flow simulations, 
e.g., that allow the tracing of a jam, need a more detailed 
description of the jamming mechanisms. For the correct 
reproduction of the upstream propagation of the down- 
stream front of a jam it is necessary to reduce the outflow 
from a jam and thus to facilitate metastable states. Here, 
the VDR model, the TOCA model and the BL model al- 
low the existence of states with a flow considerably larger 
than the outflow from a jam. 

Differences between the models can be observed in the 
jam dynamics. While the road in the VDR model is 
separated into a region with free flow and a compact jam 
that propagates upstream, the peculiarities of the update 
rules of the TOCA model lead to a jam that covers the 
whole system. 

Large compact jams appear also in the BL model since 
the slow-to-start rule of the VDR-model is included. 
However, brake lights are responsible for the generation 
of synchronized regions, i.e., regions of vehicles that are 
moving with a small velocity but high flow. 

This difference in the vehicle dynamics becomes most 
obvious in the analysis of locally measured single- vehicle 
data. On a microscopic level of description the main dif- 
ficulty lies in the reproduction of small time-headways 
that can be found at low densities in free flow and in the 
density dependence of the velocity-distance relationship. 
This important behavior of the OV-curve demonstrates 
that the driving strategy of a vehicle depends strongly on 
the traffic state while the vehicles in most modeling ap- 
proaches adjust their velocity directly according to their 
headway only, and therefore by the density. 

As a first step towards a realistic modeling of high- 
way traffic the interaction horizon of the original NaSch 
model has be to enhanced like in the TOCA, the ER and 
the HS model. However, in the TOCA and the HS model 
the cell length is not decreased which is necessary in or- 
der to reproduce realistic acceleration values. Therefore, 
the benefits of the increased interaction horizon do not 
become visible. Moreover, vehicles do react in a static 
manner to a stimulus within the horizon. In particular, 
the velocity gap matrix used in the ER model just leads 
to speed limits for certain densities. 
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A further step is the incorporation of the idea of event- 
driven anticipation. In contrast to the static reaction 
described in the previous paragraph it allows for a dy- 
namical response that will enable the vehicles to adjust 
their velocity to the actual traffic situation regardless of 
the traffic density in front. This idea is realized in the BL 
and KKW models. It turns out, that in case of the BL 
model only the introduction of brake lights, which allow 
the timely adjustment of the velocity to the downstream 
speed and can propagate in upstream direction, allows 
the reproduction of synchronized traffic. Of course, there 
might alternative ways to model synchronized traffic, but 
we believe that long-ranged event-driven interactions be- 
tween the vehicles are essential. 

The use of an effective gap by means of velocity an- 
ticipation reduces velocity fluctuations in free flow and 
leads to platoons of vehicles driving bumper-to-bumpcr. 
It is also worth mentioning that this effect is of special 
importance in multi-lane traffic as shown in [47]. 

We have seen that the BL model allows to overcome 
the problems in the reproduction of synchronized traffic 
encountered in the other modeling approaches. It re- 
produces qualitatively the observed behavior. Even the 
quantitative agreement is in most cases very good al- 
though the test scenario has neglected effects like disorder 
(different vehicle and driver types) and boundary condi- 
tions. In contrast, in most other approaches the discrep- 
ancies between empirics and model behavior can already 
be seen on a qualitative level. In particular, in the simula- 
tions of the BL model three qualitatively different micro- 
scopic traffic states are observed in accordance with the 
empirical results. The deviations of the simulation results 
are mainly due to simple discretization artifacts which do 
not reduce the reliability of the simulation results. We 
also want to stress the fact that the agreement is on a mi- 
croscopic level. This improved realism of the BL model 
leads to a larger complexity of the approach compared to 
other models of this type. Nevertheless, due to the dis- 
creteness and the local car-car interactions, very efficient 
implementations should still be possible. Moreover, the 
adjustable parameter of the model can be directly related 
to empirical quantities. The detailed description of the 
microscopic dynamics will also lead to a better agreement 
of simulations with respect to empirical data for macro- 
scopic quantities, e.g., jam-size distributions. Therefore 
we believe that this approach should allow for more re- 
alistic micro-simulations of highway networks. 

Our results show that the CA models for highway traf- 
fic have reached a very high degree of realism. The most 
complete description of the empirical findings is by means 
of the BL model. This is not surprising since the model 
has been designed in order to reproduce data of local 
measurements. But anyhow it is important to know 
which aspects of real traffic are described by a certain 
model, because in the end the aspired accordance of a 
model with empirical observations strongly depends on 
the goal of the particular application. So it is useful to 
use oversimplified model approaches in order to concen- 



trate on particular aspects of traffic flow phenomena [60] . 

Finally we want to emphasize that the results obtained 
from modelling approaches also help to improve our un- 
derstanding of the general principles of traffic flow. We 
have seen the the complexity of human behavior becomes 
more important if one wants to reproduce its properties 
more accurately. In the simplest case only the accident- 
avoidence is sufficient to reproduce the basic properties, 
like free-flow and jammed phases. For synchronized traf- 
fic, however, this is not sufficient. Here the results in- 
dicate that the dependence of the driving-strategy on 
the traffic state becomes essential. Drivers do not only 
want to avoid crashes, but also drive comfortably, e.g. by 
avoiding unnecessary large acceleration or deceleration. 
This has been emphasized in [20] and is implemented in 
slightly different form in the BL and KKW models. 

The next step would be the inclusion of other modelling 
approaches, not only cellular automata models. Using 
a different test scenario this has recently been done by 
Brockfeld and Wagner [70] . They have compared travel- 
time for various models (e.g. NaSch, VDR and OVM) 
with empirical data. Using methods from optimization 
theory to determine the best parameters it was surpris- 
ingly found that all models produce similar results that 
are not in good agreement with the data. The reason 
for this is not understood. However, the performance of 
more sophisticated models (like BL and KKW) has not 
been investigated in [70] . 
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APPENDIX A: CONTINUOUS LIMIT OF THE 
NASCH MODEL 

The adjustment of the acceleration of the vehicles in 
the original NaSch model to empirical values (that are 
about 1 m/s 2 [50]) requires the decrement of the length 
of a cell (see also [13,61,62]). 
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FIG. 34. NaSch model with different « max and homoge- 
neous initialization. 



This, however, entails an increment of the maximum 
possible velocity for a given fixed absolute value of w max 
(about 100 km/h throughout this paper). 
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FIG. 35. Fundamental diagram of the NaSch model with 
Pdec = 0.5 for different discretizations and homogeneous ini- 
tialization. 



It turns out, that already the increment of the number 
of states a vehicle is allowed to adopt leads to hysteresis 
effects of the flow. In particular, the flow can be enhanced 
in a certain density regime by initializing homogeneously 
the vehicles on the lattice compared to a pure random 
initial setup. As a result, in the limit w max — > oo [61-64] 
the system exhibits metastable states (Fig. 34) with a 



flow increasing proportional to u max , but with a rapidly 
decreasing lifetime. 

Unfortunately, increasing only w max leads to a signifi- 
cant decrement of the density of maximum flow. Thus, 
in order to keep the maximum velocity fixed, the limit 
ti mas — > oo with v max /l — const, with the length I of 
a cell has to be considered. Fig. 35 shows fundamen- 
tal diagrams for different finer discretization. Since the 
acceleration step of a vehicle is decreased considerably, 
velocity fluctuations and vehicle interactions in free flow 
are reduced. A random initialization of the system does 
not allow the high flow states so that hysteresis can be ob- 
served. On one hand, with increasing deceleration prob- 
ability pdec the stability of the homogeneous flow branch 
of the fundamental diagram decreases, but on the other 
hand the capacity drop increases. 

Unlike in the VDR model, the origin of the high flow 
states cannot be traced back to a reduction of the outflow 
from a jam but to the stability of the free flow state. 

A system with length / and deceleration probability 
Pdcc behaves like a NaSch model with cell length 1 and 
a considerably smaller deceleration probability of about 
Pdec/l- In contrast, in the congested regime the influence 
of the cell length can be neglected and a system with 
decreased cell length behaves analogously to the NaSch 
model with the same deceleration probability, e.g., the 
dynamics of the vehicles in the congested regime of the 
NaSch model is maintained (unlike in the cruise control 
limit of the NaSch model [65] where cars that are driving 
with v max have a deceleration probability Pdcc(«max) = 
0). 
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FIG. 36. Order-parameter rj for the NaSch model with a 
cell length of 1.5 m and different system sizes L for pdcc = 0.5 
and a homogeneous initialization. The inset zooms into the 
transition region. 

For realistic traffic simulations it is important that the 
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high flow states are metastable for finite systems in the 
sense that the probability for a perturbation that leads 
to a collapse of the flow is only very small. Nevertheless, 
in the thermodynamic limit the high flow states become 
unstable so that the homogeneous branch of the funda- 
mental diagram vanishes. 

In order to study the phase transition we introduced 
an order parameter that exhibits a qualitatively different 
behavior within the two phases. Because of the mass 
conservation in the NaSch model with periodic boundary 
conditions we observed the density r\ of jammed cars: 



I N 



(Al) 



In the NaSch model r) decays exponentially in the vicinity 
of the transition [66] whereas a sharp drop occurs in the 
VDR model [67] . Due to the finite braking probability in 
the NaSch model cars with zero velocity do exist even at 
densities below the transition density. In contrast, due to 
the small deceleration probability p^cc in the VDR model 
one macroscopic jam forms only at densities above the 
transition density. With increasing pd ec the transition 
smears out. 

Analogously to the VDR model, the order parameter 
of the NaSch model with a finer discretization (Fig. 36) 
shows a transition from zero to a linear dependence of 
the density. With increasing system size the high-flow 
states become unstable and the jump in the order pa- 
rameter vanishes which demonstrates the metastability 
of the high flow states. 



inequality (B2) is always satisfied if V op t(d n ) < d n — 1. 
For general A, however, this is not the case. 

This can easily be verified by initializing the system in 
a compact jam. In our simulations jams always occurred 
for global densities larger than 20 veh/km when the first 
car arrived at the jam. This simulation result has to be 
discussed in the context of the empirical results of the 
jam dynamics. Empirically one observes quite often a 
jam surrounded by free flow traffic. This includes the 
fact that cars approach the upstream front of jams with 
a rather large velocity. Unfortunately for the HS model 
these kind of configurations lead to accidents, which is in 
sharp contrast to the real situation. 

But how does one have to choose A for a given OV 
function? Using the inequalities x > \x\ > x — 1 (for x < 
0) one can derive sufficient conditions on the sensitivity 
parameter A for the model to be realistic in the sense that 
no collisions occur 



••• ma.N { r~ — 7-L — — : v > U opt (rf)| 



V opt (d) - v 
and vehicles do not move backwards 



A < min 



v opt (d) 



■ v > V opt (d) 



}• 



(B3) 



(B4) 



We checked these two conditions for the OV-function 
given in [18]. It turns out that for the chosen V op t- 
function A = 1 is the only possible choice. The up- 
per limit for A holds for a quite general class of OV- 
functions, i.e., it is the upper limit for all OV-functions 
having U opt = for some value of the gap. 



APPENDIX B: ACCIDENTS IN THE HS-MODEL: 
A STATIC CRITERION 



In order to ensure collision-free motion in a model with 
parallel update the condition 



v„(t + l) < d n - 1 + v n+1 (t + 1) 



(Bl) 



must always be fulfilled, i.e., the new velocity v n (t + 1) 
of a car has to be smaller than the number d n — 1 of 
empty cells in front plus the number v n +i{t + 1) of cells 
the preceding car moves in the next time-step. Eqn. (Bl) 
has to be complemented by the inequality < v n (t + 1) 
that ensures that vehicles do not move backwards. 

Consider now the case where the vehicle approaches 
the end of a jam, i.e., the preceding car is standing and 
will not move in the next time-step (v n+ i(t + 1) =0). 
Using the acceleration rule eqn. (4), condition eqn. (Bl) 
can be rewritten as 

v n (t + 1) + |A[V opt (d„) - v n (t + 1)]J < d n - 1. (B2) 

In order to be intrinsically free of collisions, condition 
(B2) has to be satisfied for all d and all v. For A = 1 the 
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